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General solutions for braneworld dynamics coupled with the bulk Einstein equation are derived 
under the Schwarzschild ansatz. They relate the brane metric to the exterior configurations, and 
establish fine tuning conditions for the braneworlds to reproduce successes of the Einstein gravity. 
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Einstein gravity is successful in (i) explaining the New- 
ton's law of universal gravitation for moderate distances, 
and in (ii) predicting the observed light deflections and 
planetary perihelion precessions due to stellar gravity. 
They are derived via the Schwarzschild solution of the 
Einstein equation based on the ansatze, (a) staticity, 
(b) spherical symmetry, (c) asymptotic flatness, and (d) 
emptiness except for the core of the sphere. 

There are many plausible reasons that tempt us to take 
our 3+1 dimensional curved spacetime as a "braneworld" 
embedded in higher dimensions [l]-[35]. The braneworld 
theories would be expected to inherit the successes of the 
Einstein gravity. It is, however, not true because there 
exists no Einstein equation of the braneworld itself. The 
brane metric alone is insufficient to specify the brane con- 
figuration in the bulk. Instead, we should use the brane 
position variable, which obeys its own dynamics, but not 
the Einstein equation [3]. To reproduce the Einstein 
gravity on the brane, we need finely to tune exterior con- 
figurations in general. A notable exception is the "brane 
induced gravity" [4], [12], [20], [21], [23], [30], [31], where 
the brane Einstein equation is supplied by quantum "in- 
duced gravity" [36] with a composite metric. Otherwise, 
however, closer examinations are desiered. Any entity 
can be physical only through its dynamics. Hence, we 
derive general solutions for braneworld dynamics coupled 
with the bulk Einstein equation under the Schwarzschild 
ansatz (a)-(d), and establish (fine tuning) conditions to 
reproduce the successful results of the Einstein gravity. 
We first illustrate them in (A) Nambu-Goto brane in 
asymptotically flat bulk, and extend it to (B) warped 
bulk, and to (C) Einstein-Hilbert brane. 

Let us consider a 3+1 dimensional braneworld embed- 
ded in 4+1 dimensions. Let X 1 be the bulk coordinate, 
and gu{X K ) be the bulk metric at the point X K [37]. 
Let the brane be located at X 1 — Y I (x^) in the bulk, 
where x M are parameters. The dynamical variables of 
the system are Y I (x^), gu(X K ) and matter fields, but 
not the brane metric g^ u = Y I Y'lgij(Y K ) [37]. Here we 
consider the simplest models given by the action, 
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RP »u\ and rL uk are the 



-g{\ + 1 R)d i x~ U^§(X + jR)d 5 X + S m , (1) 



brane and bulk curvature tensors written in terms of g^ v 
and gjj, respectively, S m is the matter action, and A, 7, 
A, and 7 are constants. We add no artificial fine-tuning 
terms. The equations of motion derived from (1) are 



-Ay'"' + - j R» v - X -Rg» v ) - V" 



-^9u + j [ Ri.j - 2 R 9 IJ 



T, 



IJ, 



= 0, (2) 
(3) 



and those for matters, where T M „ and T/j are the en- 
ergy momentum tensors with respect to g^ v and gu, 
respectively [37]. Eq. (2) is the Nambu-Goto+Regge- 
Teitelboim equation [3] for the braneworld dynamics, and 
eq. (3) is the bulk Einstein equation. The brane Einstein 
equation does not hold. Hence, we solve (2) and (3) at 
the brane, and examine their implications. (If necessary, 
they can be continued into the bulk in known methods 
[24]-[29], [35].) We first illustrate them in a simple model 
in (A), and discuss extensions in (B) and (C). 
(A) Nambu-Goto Brane in Asymptotically Flat Bulk (7 7^ 
0, A 7^ 0, A = 7 = 0) We introduce the Gaussian normal 
coordinate (A' 1 , A 4 ) = (x 1 - 1 ,^) with £ = at the brane, 
and expand the bulk metric gu in £ as 



9nu 
9ni 



o(a, 



344 = -l + 0(£ 3 ), 



(4) 



where and C^„ are functions of x^. The former 
is the extrinsic curvature: K^ v — — rt^|j=o with 
the bulk affine connection T JK - We introduce the time, 
radial, polar and azimuth coordinates, x a = t, x 1 = r, 
x 2 = 6, and a; 3 = 4>, respectively, according to ansatze 
(a) and (b), which imply that gu is diagonal, and goo, 
922 = 933/ sin 2 8, and 544 are functions of r and £. 
We parameterize them with 



= diag(/, h, r 2 ,r 2 sin 2 



(5) 



= diag(a, b, c, c), Uj? = diag(u, u, w, w), (6) 

where Ujf = (R 41 * v4—o~v R AX m)\£=o represents the degrees 
of freedom of C^ u = K*K\ V + R* ^4^=0, an d /, h, a, b, 
c, u, v, and w are functions of r. Then, eq. (2) implies 



where g — detg^, g — detgu, R — R^ , R = Wj, 



a + 6 + 2c = 0, 



(7) 
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while eq. (3) at the brane leads to 

-(h-iy/r + (1 - h-^/r 2 - 26c - c 2 + u = 0, (8) 

-f'/rfh + (1 - h-^/r 2 - 2ac - c 2 + v = 0, (9) 
-((A)- 1/2 /')'/2(A) 1/2 - (//M'/2rf 

— a6 — ac — be + w = 0, (10) 

-mr 1/2 f'YWh) 1/2 - (f/h)'M 

+ (1 - /i _1 )/r 2 -ab- 2ac - 2fec - c 2 = 0, (11) 

where ansatz (d) is used, and contributions from the 
brane action to Tjj are neglected. We have five equa- 
tions (7)-(ll) for eight variables /, h, a, b, c, u, v, and 
w. Accordingly, the general solution involves three arbi- 
trary functions, for which we here choose a, c and v. 

We first eliminate /, /', /" and b from (7), (9) and 
(11) to obtain the differential equation, 

2r 2 (4V + rV -A) 



(rh- l y = 1 + 



r 2 V 



h(l + r 2 V)-3 



with V 



2ac 



A = 2a 2 + 4ac + 6c 2 , 
To solve eq. (12), we rewrite it into the form 

Z' = P + Q/(l + Z/r) 

with Z = 3(/i _1 - l)/2 - Vr 2 /2, 
P= {A-AV-Wr/2)r 2 , 
Q = (1 + VY 2 )(i4 - 4V - V"r)r 2 /2. 



(12) 

(13) 
(14) 



(15) 
(16) 
(17) 
(18) 



A sufficient condition for existence of the unique solution 
to (15) with the initial value Z = Zr \ at p = 1/r = is 
that, in a region D in the p-Z plane including (0, -Z(o)), 

P/p 2 + Q/p 2 (l + Zp) is continuous, and (19) 
3M : real V(p, z x ), (p, z 2 ) G D 

|Q/p(l + zip)(l + 2 2i o)| < M, (20) 

which we assume here. Then, the solution to (15) is given 
by Z = lim„_ 



Z( n ) with 



J (n) 



J (0) 



f 



Qr 



r + Z(n-1) 



dr. (n > 1) (21) 



Another convenient form of the solution would be given 
by series expansion in 1/r. Eqs. (19) and (20) indicate 
that Pr = —Qr =finitc at r = oo. Hence, we assume 
existence of (at least, asymptotic) expansions: 



Z = Z + Z 1 /r + Z 2 /r 2 + 
P = P 1 /r + P 2 /r 2 + -.., 
Q = Qi/r + Q 2 /r 2 + •••, 



(22) 
(23) 
(24) 



where Z i} Pi and Qi (i = 1,2, ■■■) are constants, and 
Zq = Z( y Then, eq. (15) implies the recurrence, 



z _ _ y-y Z k [(n - l)Z n - k -i + 2P w -fc] 



J n+1 



fc=0 



2n 



(25) 



for n > 1, where S n = P„ + Q n , or, explicitly, 

Zi = -Z P 1 - S 2 , (26) 
Z 2 = [Z 2 Px - Z (P 2 - Q 2 ) + PiS 2 - E 3 ]/2, (27) 
Z 3 = -[2Z 3 P! + Z 2 (5P 2 + 2Q 2 ) 

+Z {SPl + AP^ + 1P X Q 2 - 2Q 3 ) 
+Pi 2 E 2 -PiE3 + 2Q2S2 + 2S 4 ]/6, (28) 



Once the function Z of r is thus determined, the general 
solution of the braneworld is given by 



h = l/(l + 2Z/3r + Vr 2 /3), 
h(l + Vr 2 )-l 



f = exp 



/ 



dr, 



u = [(r/h)' - l]/r 2 - 2ac - 3c 2 , 
w = —(u + v)/2, 
b = —a — 2c, 



(29) 

(30) 

(31) 
(32) 
(33) 



with three arbitrary functions a, c and v. The braneworld 
dynamics (2) is crucially important. If we omitted it, and 
used only the bulk Einstein equation (3), we would have 
had the wrong solution without (33) and with 

A = -2ab - Aac - Abe - 2c 2 , (34) 
u = [(r/h)' -l]/r 2 + 2bc + c 2 , (35) 

in the places of (14) and (31), respectively, which would 
cause large differences in the consequences below. 

Now we compare our general solution with the results 
(i) and (ii) of the Einstein gravity in 3+1 dimensions. 
The latter is given by (8) and (9) with a = b = c = u = 
v = 0, and is solved to give the Schwarzschild solution 
/ = h^ 1 = 1 — ji/r with an arbitrary constant p, which 
is determined phenomenologically. In geodesic motions 
of objects, the function Vn = / — 1 = — p/r serves as the 
Newtonian potential for gravitational forces for moderate 
distances, as renders a precise explanation for the New- 
ton's law of the universal gravitation. The angles A</?e 
and A^e of light deflection and planetary perihelion pre- 
cession, respectively, due to the stellar gravity are given 
by A(^e = 2/x/ro and A^e = 3np/(l — e 2 )a, where ro 
is the minimum of r on the light orbit, and e and a are, 
respectively, the eccentricity and the average distance of 
the planetary orbit. These predictions precisely match 
with the observations. 

Now we turn to our general solution. The Newtonian 
potential Vn = / — 1 has large arbitrariness owing to 
those of a, c and v. In fact, we can always choose ap- 
propriate a, b, c, u, w and v for arbitrary / and h with 
R > 0, while we have no solution for R < 0. Therefore, 
we cannot predict that Vn = —p/r, unlike the Einstein 
gravity. The singularity structures also depend on exte- 
rior configurations. The conditions (19) and (20) would 
be useful in investigating them. 
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In ordinary observations with r » ji, the fi/r term 
dominates Vn • Then, we parametrize the deviations from 
the results of the Einstein gravity by 



/ = 1 - fi/r + h/r 2 + h/r 3 + 
hr 1 = 1- n/r + e 1 /r + £ 2 /r 2 + ■■ 



(36) 
(37) 



with constants ft (i > 2) and e, (i > 1). The asym- 
totic flatness implies f,h—>l and /', h! — > as r — > oo. 
Then, eq. (9) (with A = 0) implies Vr 2 0. This and 
cqs. (17)-(20) indicate Ar 2 — > 0, and therefore (34) imp- 
ies ar, br, cr — > 0. These and Vr 2 — > indicate vr 2 — > 0. 
Hence, we assume existence of (at least, asymptotic) ex- 
pansions of the arbitrary functions: 



a = a 2 /r 2 + a 3 /r 3 + ■ 
c = c 2 /r 2 + c 3 /r 3 + •• 
v = v 3 /r 3 + v^/r 4 + ■ ■ 



(38) 
(39) 
(40) 



with constants a^, Cj (i > 2) and Vi (i > 3). Then, the 
present general solution indicates 



£l = V3, 

f 2 = 3/iw 3 /4 - Vi - a 2 2 
£ 2 = fiv 3 /2 - Vi - 2a 2 2 



2c 2 2 , 

2a 2 c 2 - 5c 2 2 



(41) 
(42) 
(43) 



In terms of them, the angles Aip and Aip of light de- 
flection and planetary perihelion precession, respectively, 
due to the stellar gravity are given by, up to 0(r 3 ), 

Acp « A^ e (1-£ 1 /2 M ) = A<^e(1-^3/2 / u), (44) 
AV- « A^ E [l-( M e 1 + 2/ 2 )/3 M 2 ] 

2c 2 2 )- 



= A^E 



1 _5v3 + 2(f 4 
6fi 



a 2 



3^ 2 



(45) 



Consequently, the condition to match the observations is 



v 3 pa 0, Vi + a 2 + 2c 2 



0. 



(46) 



The condition for the bulk Einstein equation (8-11) at 
the brane to imply the brane-Einstein equation would be 

u = 2bc + c 2 , u = 2ac + c 2 , w = ab + ac + bc. (47) 

They are, however, not all physical. Under the dynamics, 
(47) leads to a 2 +2ac+3c 2 = 0. Hence, the brane Einstein 
equation holds only when 



w 



0. 



(48) 



(B) Warped Bulk (A ^ 0, 7 ^ 0, A ^ 0, 7 = 0) The 
dynamics (eqs. (2) and (3)) has a solution with a flat 
brane at |£| < S in empty warped bulk [19]: 

g u = diag(F, —F, -Fr 2 ,-Fr 2 sin 2 9, -1), (49) 

where S is an infinitesimal constant, and F is a smooth 
function of £ such that F'| ?=0 = and F = e~ 2k ^\ 



for |£| > 6 with fc 2 = —A/67. We seek for the gen- 
eral solutions, where the metric gu approaches (49) as 
r -)• 00. Hence, f^| ?= ±«5 -> Tkg^, f^| ? =o -> 0, and 
i? 4M !y4|^=±5 — > fc 2 f>£J. The ansatz (c) is retained on the 
brane, but (d) is slightly violated by matter distributed 
in |£| < S. We assume that the brane-generating interac- 
tions are much stronger than gravity at short distances 
of 0(5), while their gravitations are much weaker than 
those by the core of the sphere. Then, T M „ in |£| < S be- 
comes independent of r, and so does (^^=±5 — r^J{=o) 
[38], so that it is equal to its asymptotic value ^fkg^. We 
define the functions a, b, c, u , and w 1 * 1 of r by 



r^|f=o = r^|| =±5 ± kg ^ = diag(a, b, c, c), (50) 

Xi)k=±S 

= -3k 2 S£ + diag(u ± ,5 ± ,w ± > ± ). (51) 



(R 4 ^i-S^R 4X X i)\^ ±s 



Then, eq. (2) with 7 = and eq. (3) at £ = ±5 imply the 
same equations as (7)— (11) in (A) with a, b, c, =F 2ka, 
=p 2kb and =F 2fcc in the places of a, b, c, u, v and 
w, respectively. Therefore, we obtain the same results as 
in (A) above, though the geometrical interpretations are 
different. 

(C) Einstein Hilbert Brane (7 ^ 0, 7 ^ 0, A = A = 0) 
The brane dynamics (2) implies, instead of (7) in (A), 

[-(h-iy/r+il-h- 1 )/^ 
+[-f'/rfh+(l-h- 1 )/r 2 ]b 

+2[-((fh)-^ 2 f'Y/2(fh) 1 / 2 (f/h)'/2rf}c = 0, (52) 

while cqs. (8)— (11) are unchanged. If we eliminate /, /', 
/", b, u, and w from eqs. (8)— (11) and (52), we obtain 



(rh- 1 )' = l + r 2 V 

r 2 [rV'(V - 2ac - 4c 2 ) + 5V 2 - 2y a (2a + c)] 
+ h(l + r 2 V)(V - 2ac - 4c 2 ) + 2(a + 2c)(2o + c) ' 



(53) 



Eq. (53) has the same structure as (12) as far as h is con- 
cerned. Hence, we can apply the same method as that 
in (A) above. In the present case (A = 0), the condition 
(47) automatically implicates the brane dynamics (52). 
Therefore, (47) just gives the condition for the brane Ein- 
stein equation to hold, but not (48). 

In any cases, the successful results of the Einstein grav- 
ity are reproduced through fine tuning, and are not "pre- 
dicted" by the braneworld theories. It would be an urgent 
and important challenge to seek for physical foundations 
of the fine-tuning conditions. 
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